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ABSTRACT 


Application of radiation and first order chemical reaction to the behavior of fluids in electromagnetic fields is investigated 
analytically. The study reveals that the behavior of fluids in electromagnetic fields is governed by electrical conductivity which gives 
rise to electromagnetic forces and mechanical forces. The effect of electrical conductivity decreases the magnetic field and increases 
the electric field while increase in the electrical resistivity, increase magnetic field strength and decrease the electric field strength. 
Critical observations also reveal that increase in the flow of energy, decreases the heat transfer coefficient and increase in the 
chemical reaction term results in a decrease in the mass transfer coefficient. 
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1. INTRODUCTION 


Fluid in physics is made up of liquids and gases. Life as we know it would not exist without fluids and without the behaviour of fluids. 


Essentially, when fluid is introduced into an electromagnetic field, it is primarily governed by electrical conductivity. This effect gives 
rise to not only electromagnetic forces but mechanical forces as well. An alternating magnetic field applied to a conducting fluid will 
induce electric currents in the fluid and hence a Lorentz force distribution. The Lorentz force is rotational and is set into motion. This 
assertion was corroborated by Branover [1] and he went further to say that electromagnetic field exists everywhere on the earth's 
surface and several studies on its influence are abounded. In a work Published by Eisenschitz and Cole [2], they opined that if an 
insulating liquid is placed between the plates of a condenser, the electric field gives rise to induced intermolecular forces which are 
proportional to the square of the field strength. Symons [3] investigated that fluid is unique but not anomalous and reported that, it 
appears that electric and magnetic fields have opposite effects on fluid clustering. Colic and Morse [4], in the elusive mechanism of 
the magnetic memory of water, suggested that, in addition to the breakage of hydrogen bonds, electromagnetic fields may perturb 
the gas/liquid interface and produce reaction oxygen species. Zhang et al [5], studied the behaviour of dielectric liquids in electric 
field and reported that, a silicon oil is subjected to electric field and the result using the Ohnesorge number revealed that due to the 
low electrical conductivity and low electrical relaxation time of the silicon oils, only unsteady transient jets were found. Ninno and 
Castellano [6], in their findings, reported that extremely low frequency electromagnetic fields have significant and lasting effects on 
liquid water. Using a weak field, adjusted to give a magnetic field of 45 micro Tesla on glutamic acid solution cause changes in the 
PH shifting towards the de-pronated species. Sherwood [7] in his study of fluid droplets in electric and magnetic fields reported that, 
a fluid held by surface tension will deform when an electric or magnetic field is applied and depend largely on 
permittivity/permeability and conductivity of the drop and of the surrounding fluid. Zanh and Pioch [8], reported that analyses and 
measurement have shown anomalous behaviour of ferro fluids in alternating magnetic fields, this they investigated by pumping 
ferro fluid which the flow direction can reverse as a function of magnetic field amplitude, frequency and direction. In a research 
carried out by Daniel [9], a new generation of magnetic sensitive ionic liquids-magnetic ionic liquids, exhibiting a strong response to 
a magnetic field was developed. This shows that its physical properties are influenced by the presence of an applied magnetic field. 
Several studies are also abounded in the behaviuor of fluids on non electromagnetic forces. Few are Prasad and Bhaskar [10], 
Ferdows et al., [11], Hayat et al., [12], Subhakar and Gangadhar [13], Shivaiah and Anand [14] and Gnaneswara and Bhaskar [15], all 
studied the effect of similar and different non-electromagnetic forces on the behaviour of different types of fluid and made far 
reaching suggestions. The aim of this study is to consider the effect of electromagnetic forces and mechanical forces on fluid flow 
with emphasis on radiation and chemical reaction as non electromagnetic forces. This combination of forces caused by electrical 
conductivity on the behaviour of fluid in electromagnetic field is an extension of the literatures cited. 


2. FORMALISM 


The equation of continuity 


op 
——+¥. =0 (1) 
ar 1") 
The momentum equation 
0 
pa =—-Vp-— (JXB)+ LN’?v + pg (2) 
t 
Maxwell’s equations 
B 
VxE = = (3) 
Ot 
VxH =J+ 2 (4) 
Ot 
Ohm's law for a one component conducting fluid 
J = o(E +vxB) (5) 
Experimental law of Gauss for magnetic field 
V.B=0 (6) 


Supplementary equations (constitutive relations) in isotropic medium 
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D=6&,E (7) 
1 

H=—B (8) 
Ho 


The energy equation 

oC oT ,OT %, 
~ 2 

"dy dy” Oy 

The concentration equation 


2 
pC, = DE -kC (10) 


Considering the Boussinesq approximation in the body force to incorporate the temperature dependent density and concentration 


(9) 


dependent density, equation (2) can be written as 
Ou . 
p=, =—Vp ~(JXB)+ ui + pg BT -T))+ pgB*(C-Cy) (11) 


with the boundary conditions 


u=1,T =1,C=1at p=0 (12) 
u>0,T>0,C 30 at pol 


In order to consider the effect of radiation on an optically thick model in which the thermal layer becomes very thick or highly 
absorbing as described by Rosseland approximation, Cogley et al [16] as 
é 4¢ O°T* 

dq, __ 4G 4a 


oy = 3a Gy” 


where (,, is free charge density, J is free current density, E is electric strength, B is magnetic induction, H is magnetic strength, &, is 


r 


permittivity of free space, Ly is permeability of free space, u and V are velocities and Xx ,¥y directions, t' is time, C is fluid 
concentration, p is pressure, 9 is density of fluid, 2 is viscosity of fluid, O is electrical conductivity of fluid, B is imposed 
magnetic induction, g is acceleration due to gravity, @ is thermal expansion due to temperature, Bg is thermal expansion due to 
concentration, T is temperature of fluid, T,, is free stream temperature, C,, is free stream concentration, k is thermal conductivity 
of fluid, C, is specific heat at constant pressure, g,. is radiation term, C. is concentration susceptivity, ie is mean nanofluid 


temperature k? is chemical reaction term, D is chemical molecular diffusivity, ¢ is the Stefan-Boltzmann constant and @ is the 


absorption coefficient. If temperature difference within the flow of the fluid is sufficiently small, we can approximate T* using 
Taylor series expansion about 0 and obtain 


T* =4T3T - 37, (14) 

Using equation (3), the time dependence of magnetic field is 

OB 

— =-(VxE (15) 
ap 7 NE) 

Equation (5) is put into equation (14) and the result is 

OB J 

—_ = —V(vxB)- v=] (16) 
Ot o 


Neglecting the displacement current in equation (4) and using equation (7), equation (15) becomes 


oF V(vxB)+ | ve(vxB) (17) 
Ot ie 
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Simplifying, using equation (6), equation (17) reduced to the form 


a ~ V(vxB)+—_V?B (18) 
Ot Myo 


If the fluid is at rest, then equation (18) reduced to 


OB __1 vp (19) 
Ot MWe 


with the boundary conditions 


B(0,t)= B(z,t)=0 B(y,0)=1 (20) 


Method of solution 
One dimension of equation (19) is assumed and taking the finite Fourier Sine Transform on both sides, the expression takes the form 


a n2 
e BUy.1) Sinnydy (21) 


oy 
where n is an integer 


Let By (n,t) = [ Bly, 1)Sinnydy, then 
0 


dB,(n,t) ¢OB(x,t) ... 1 7O°B(x,f) .. 
a(t) =| ( eer | ( ) sinnydy (22) 
dt ‘ Myo 5 
Simplifying and imposing the boundary conditions, equation (22) reduced to 
1-—Cosnx 2 
B,(n,t) = [ee nm ‘ (23) 
Hyon 
Applying the inversion formula for finite Fourier Sine Transform, following Gupta [17], equation (23) is written as 
2 S| { 1-—Consxz afte Nhe 
B(y,t) = >» Joos " |Sinny (24) 
n=l Hyon 


If the conductivity is very large or infinite, equation (18) reduced to 

OB 

— = V(vxB) (25) 
Ot 


To determine an expression for the electric field intensity, equation (5) is put into equation (15) and the expression takes the form 


1 10E 
1a) 1 yp - 


ov Ot v Ot 


a is taken on equation (4) and also taking into consideration VxB = Mot as well as the constitutive relation of equation (7), 
t 


equation (26) will transform into 


OE ovoE 
an areal (27) 

Ot E, Oy 
with the boundary conditions E(y,0) = E(y, 1) =] E(0,1) =] (28) 


Using the method of separating the variables and imposing the boundary conditions of equation (28), the solution of equation (27) 
is written as 
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E(y,t) = Cosk ,tExp — (22429) (29) 
ov 


where k, is a constant 


Equation (25) shows that the magnetic lines of force are in motion with fluid velocity and that the magnetic flux through any circuit 
moving with fluid velocity is constant in time. Also, equation (11) and the use of equation (5) can be written as 
Ov 
ar =F- o(E + (vxB))xB (30) 
t 


where F is the sum of all the non-electromagnetic forces only. Simplification of equation (30) and for fluids of infinite conductivity, it 
reduced to 
Ov 


—=F (31) 
Ot 


This shows that flow parallel to B is governed by non-electromagnetic forces only, which is the mechanical behaviour of the system. 


To tackle the effect of mechanical forces imposed on the system, equations (9)-(12) are transform into dimensionless form. The 
dimensionless quantities used are 


T-T C-C 
(j= 2.9= = tReet oe = C= 7 
xt Xx uy VoP 1. ed C, =O 
y k? 16¢<T; 
pra Cr)  _ p sea g = N= os (32) 
k Ox D v 3a C, 
eB(T -T,)u eB (Cc, )u 
Gr, = ——.__. Gr, = ——_—_ 
Vout Vou 


By using the dimensionless quantities, the governing hydrodynamic equations are transformed into the form 


Ou 1 Ou 
one ee ay? + Gr,8 + GroC (33) 
e 
2 
pa peat NO (34) 
r 
2 
a == ao Kae (35) 
Cc 


and the boundary conditions are 


u=1,0=1,C=l1at y=0 (36) 
u>0,0>0,C>0 at yol 


where Re is the Reynolds number, Pr is Prandtl number, Sc is Schmidt number, Gr, is Grashof number due to temperature, 
Gr is Grashof number due to concentration, N is dimensionless radiation term, ky is dimensionless chemical reaction term, @ is 


dimensionless temperature, C is dimensionless concentration. 


Solution Technique 
The solution to equations (33) — (35), assumed the form 


u=U(y)e", 0=P(yje™, C=w(yje", P=P(yje™ (37) 
Substituting equation (37) into equations (33) - (35) reduced to 
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Re 'U"(y) +nU(y) + Gr, @(y) + Grew(y) — P, =0 (38) 
Pr! @"(y)+(n—-N)d(y) = 0 (39) 
Sc'y"(y)+(n—-k, wy) =0 (40) 


with the boundary conditions 
U =e”, g=e", y =e” at y -0 
U->0,¢>0,w>0at yol (41) 


To solve equations (38) - (40), Homotopy Pertubation Method as used by He ([18] — [20]) is applied as follows 
HU, p)=(1-p)Re'U"(y)|+ p[Re U"(y) + nU(y) + Gr, @(y) + Grew(y) - P, |= 0. (42) 


H(g, p) =(1— p)|Pr! ¢"(y) + p[Pr! 6"(y) + (n— Noy) |= 0 (43) 
Hy, p) =(1- p)|Se"w"(y) + plSctw"(y) + (nk, wy) |= 0 (44) 


Assume the solution to equations (42) -(44) be written as 


U=U,+ pu,+pU,+ ates Satetenes 
— = by + PO, + P'Dy Hevversvee ss 
Vo =Wot DW t+ DW terereee 


Substituting equation (45) into equations (42) — (44) and comparing the coefficients p’, p', Pp and Pp , the results are 


p’: Re'U}-P.=0 (46) 
p’: Pr'¢=0 (47) 
p’: Sc'wi=0 (48) 
p': Rel Us+nU, +Grrdy) + Grp =0 (49) 
p': Pr'¢’+(n—N), =0 (50) 
p': Sc'wt+(n-k, yw, =0 (51) 
p: Re'U%+nU,+Gr,¢,+Gr,y, =0 (52) 
p?: Pr'gf+(n—N)d, =0 (53) 
p’: Sc'wht+(n-k, Ww, =0 (54) 
p?: Re'U%+nU,+Gr,¢,+Gr,y, =0 (55) 
p: Pr' @%+(n-N)p, =0 (56) 
p?: Sc'wl+(n-k, w, =0 (57) 


with the modified boundary conditions 


U, 20, 30,3. =2 a 7=0 
% =9%, =>, =, e”" at y=0 
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at 


at y=0 (58) 


Wo =Wi HW, =W3 =e 
U,=0, =0, 0,0 a. yeu" 
$, =% =, =, =0 at y=e™ 

Wo HV =V2 =W, =0 at y=e™ 


Solving equations (46) — (48) and imposing the appropriate boundary conditions from equation (58), the results are 


1 1 
U, =-—ReP.y* +—ReP.e” y-1+e™ (59) 
2 2 
@) =-yte™ (60) 
Wo =-y+e™ (61) 


Similarly, the solution to equations (49) — (51) after imposing the appropriate boundary conditions from equation (58) are 


1 1 1 
U,= ~yRenUoy’ ~ 5 ReGrgyy” ~ 5 ReGrevoy" 


1 
+5[RenU ye" +ReGr,g,e" +ReGr.w,e™ ~2}y +e” (62) 
1 2: 1 nt nt 
b =—>Pr{n— Nyy + 5 Pela Nive —lly+e (63) 
1 2 1 nt nt 
Y, =~ 5 Seln— ko Woy + 5 5cln—ky Woe —llyt+e (64) 


Further, the solution of equations (52) — (54) and the use of the appropriate boundary conditions from equation (58) are respectively 


1 il 1 
U,= ~zRenuy” ~ 5 ReGrgy" ~gReGrewiy" 


+ (RenU,e" + ReGr,g,e" +ReGroy,e" ~2)y+ e” (65) 
1 2, 1 nt nt 
b, =—5Pr(n— Nay + 5 Pla Nhe -lly+e (66) 
1 2 1 nt nt 
V2 =—5Scln—ky Wy + 5 5cln— ko ie —llyt+e (67) 


Finally, following the same procedure, the solution of equations (55) — (57) are 


1 1 1 
C= ~zRenUay ~ 5 ReGrrry” ~gReGrevay" 


1 
+ - (Re nU,e”" +ReGr,¢,e" +ReGrow,e” — 2)y +e™ (68) 
1 2 1 nt nt 
@, =—5Pr(n—N)bsy > 5 Pn N)b,e all yate (69) 
1 2 1 nt nt 
VY; = 5 Seln— ko Woy + 5 5cln— ko Woe =] eke (70) 
Equations (58) — (69), give the expressions for the velocity, temperature and concentration as 
WSU Uy PU ge achieves (71) 
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PD = Py +B, +o FPyrescccreees (72) 
SY PV Fs PY osciotes (73) 


Equations (71) — (73) is put into equation (37) and the complete solution takes the form 


U = Uy $y $y + Uyeccceees (74) 
2 = 6,40, HOY Pej iscasise (75) 
Oe Os Or gl Ora rere (76) 


Skin friction, Nusselt number and Sherwood number 
The Skin friction is obtained as 


c= ey = 5RePle" y - =e" + 3(e" y + 5 ReGr, (e" ii + =Re Gr. (e" y - 5Re ne" y (77) 


The rate of heat transfer coefficient as 


n= {B) darian) +Ley Merl eM O-Nlerl ee 


and the rate of mass transfer coefficient as 


Sh= {<| =4+4Sc(n—k, Ye" ij + a(se) (n—k,) (e" y + a(se) (n—k,) (e" y’ (79) 
y J 9 


3. RESULTS AND DISCUSSION 

fy =4md0" NIA? x6, =8.85x10 °C’ Nm’ tf =1,0=1,P =1 

o(Qm)" = 2,4,6,8,10,0 | (Qm) = 2,4,6,8,10, N = 2.5,3.5,4.5,5.5,6.5 

ky =1.8,3.6,5.4,7.2,9.0, Gr, =1.5,2.5,3.5,4.5,5.5, Gro = 2.8,3.8,4.8,5.8,6.8 


Figure 1 The dependence of Magnetic field strength on Coordinate with electrical conductivity term © varying 
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1 20000 - 


1 40000 


-1 
Figure 2 The dependence of Magnetic field strength on Coordinate with electrical resistivity term 2 varying 


-1 
Figure 4 The dependence of Electric field strength on Coordinate with electrical resistivity term 2 varying 
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Figure 5 The dependence of Skin friction on Reynolds number with Grashof number due to Temperature Gr, varying 


30 | 1015 
251 1015} 
— 201 1015} 


151 1015 | 


101 1015 | 


Figure 6 The dependence of Skin friction on Reynolds number with Grashof number due to concentration Gre varying 


Nu 


Figure 7 The dependence of mean Nusselt number on Prandtl number with Radiation N term varying 
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Figure 8 The dependence of Sherwood number on Schmidt number with chemical reaction © term varying 

Figure 1 shows a graphical relationship between the magnetic field strength and the space coordinate as the electrical 
conductivity of the fluid varies. The physical explanation of the graph is that, increase in electrical conductivity of the fluid 
corresponds to a decrease in the magnetic field strength. This observation appears to be the reverse case as depicted in Figure 2, 
which shows an increase in electrical resistivity of the fluid result in an increase in magnetic field strength and later dropped to zero. 
The effect of increasing electrical conductivity as shown in Figure 3 reveals that, the electric field intensity is increased and the 
reverse is the case with Fig4 as the electrical resistivity is increased. Comparison of Figure 1 and Figure 3 as well as Figure 2 and 
Figure 4 showed that the electric and magnetic fields have opposite effects, their reciprocals also and this laid credence to the 
findings of Symon [3]. The relationship between the skin friction and the Reynolds number as depicted in Figure 5, shows that 
increase in the Grashof number due to temperature correspond to an increase in the skin friction of the induced electrically 
conducting fluid and this observation is also the same with Grashof number due to concentration as depicted in Figure 6. A close 
observation of Figure 7 reveals that an increase in the radiation parameter, results in a decrease in the transport of heat coefficient 
considering the atmospheric Prandtl number of air. The explanation of this observation is that, an additional heat source due to flow 
of energy may lower the heat transfer coefficient of the fluid. Increase in chemical reaction which led to improved concentration of 
the fluid causes a decrease in the mass transfer of the fluid owing to slow motion as a result of increased velocity as shown in Figure 
8. 


4. CONCLUSION 


When a fluid is introduced into an electromagnetic field, two major forces come into play as a result of the electrical conductivity of 
the fluid. The inclusion of radiation and chemical reaction terms to the existing non electromagnetic forces and to investigate the 
flow of the electromagnetic forces which have opposite effects are also interesting and conform to existing results. More study on 
the combined effect of electromagnetic forces and mechanical forces with the inclusion of other non electromagnetic forces of 
interest shall be our focus in the next study. 
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